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ABSTRACT 

J? SJ? muon groups observed at Baksan underground 
scintillation telescope is analysed. In this analysis 
we compare the experimental data with calculations 
J? 8ed ° n a superposition model in order to obtain ’ 

S e ,H e f feCtlVe at ° mlc n ff >er ig f P rima ry cosmic rays 
in the energy range lO^-io^ev. 

— Patroduct ion . Our general approach to the problem was nrp- 

Baksan in (If ?hf first experimentSlte” Tron™ 

Baksan in (2). The first attempt to analyse Baksan multiolici- 

galorf f f Now fT of , priraar y composition was made in Ban- 
galore (3). Now we have to correct one of conclusions made in 
(3) concerning effective atomic number. The change came for 
two ™8ons:l)now a smaller zenithangle intervfl was accep- 
• , 30, corrections for individual events inside this 

interval were made. This affected the obtained decohe?enee 
curve, making it 20* narrower, 2) the second reason is the choi- 
ce of more reliable function N u =F(E n ) which represents the 
mean number of muons from primary nucleon of energy Eq . 

2 . Experiment . The observed multiplicity or muon number spec- 

thf^dif f eren" “-J 1 ?- ^Vl Practically the same as 1^(3) 
u,^ f > dlfferen ? ’ 'Or m 3.3 being not signif icant . All the 
differenc e with (3) is in calculated curves. The experimental 
decoherence curve is shown in fig. 2. It deviates slightly from 

curve e is°used 1 to of Wlt J P^ameter (12*1) m.The defoherence 
curve is used to choose the parameters of assumed lateral 

istribution functions (LDP) in such a way as to fit the 
experimental curve, D(R) . F^ (r )-s-F/j (r) : examples of tried LDF 
Then the important function A (R) (see below) does not 
depend much on the choice of the type of LDF. 

3 _. Calculat ion s of expected muon numbe r spectra for a given 
been a made:' To do this followi ng assumption Sa ve 

1. the muon group from A is a superposition of A independent 

groups from constituent nucleons, P 

2. LDF does not depend on A and energy per nucleon En, 
^.fluctuations of the muon numbers are pure Poissonian 

4. the power law energy spectrum of primaries of a given A 

5. the mean number of muons of energy E u ^E fh from primary’ 

nucleon of energy Eq was taken as — 0.s>5s(^£f > ‘^ 7 -o.SSj) i ' ss 
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for E th = 0.2,2 Tev, Eg in Tev on the basis of numerous Monte- 
Carlo simulations of cascades in atmosphere (5,6). 

The problem is presented 
by a system of equations: 

A/=A -fee.) 

A/' = //•&( R) 


- \e. ^ / 4 

I A C m ) = 3 5 ' 


Here: N is the mean number 
of muons E >-Ex. h from nuc- 
leus A with total energy 
A«Eg;A(R) is the fraction 
of N which is expected to be 



l 5 10 2 

Fig. 1 

inside detecting area when the core of group is at the dis- 
tance R from telescope centre (A(R) is calculated using as- 
sumed LDF) ;lp A (/V' ) is distribution of N‘ (expected numbers in- 
side area), which is obtained by solving eqs./l-3/i JA( m ) i s 
the intensity of multiplicities m for a given A; I A (m) the 
same, but normalized to the total muon flux through detecting 

1 1 1 ‘ 


area. The resulting from these 
calculations X^(m) are shown 
in fig. 1 for several A and 
mixed composition . Obviously , 
these results depend on ^=1.7 
(which is in agreement with 
experiment) and on functions 
f(En) andA(R). 

4 .Discu ssion. To demonstrate 

& "P r to 

>5 


D (R) 


_ [ -0.7 _ 

F , = J 


1 - c -coo © 
r S« 

.22 Tev 



1- 10m .5 

2 - 11 .6 

3- 12 .75 

4- 6.2 


the sensitivity of A 
assumed functions f(Eg') and 
A(R) let us use "Ng-appro- 
ximation" , introduced in (1),' 
valid for m^ljwhich consists 
of following: g 

1. at large N f (Eg ) ~ Eg , 
then 

at large N’, 

2. suppose this power law distribution is valid for 
and is cut to zero for N’< N A ; determine N* by 

so providing the collect total number of muons. 

In a similar way we determine Nq which corresponds to the 
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case A(R)=l.Fig.3 shows the original f(E n ) (5) 
dence of number of muons N, u 

Eu^>220 Gev on the primary 
nucleon energy Eg and its 
power law version in "Nq- 
approximation" .The total 
number of muons after integ- 
rating over primary spect- 
rum should be the same for 
both versions .Substituting 
the functions f(Eo) and 
^k(N') in ({1/ and / 4 / by 
power law versions using the 
mentioned procedure we ob- 
tain a simple solution of 
eqs. /1-5/: 

<ZP y 

& Oil dti 


the aepen- 


„ o.A7 1.55 

|CEo)=Q355 (Eo -0 - 5 ^) ‘ 


IX"0 = 


if-J 

5 


ml 





E^^-0.22 Tev 


E n (Tev/n) 

■ q » I 




Pig. 3 


s 

<M 


- IS 


where N£=N 0 • A-Aeff > A e f f = R d R (A 0*0) ^ 

the "effective fraction" of total number of muons in the group 
expected in detection area. For m-»N^ and mas>l 

IM = (Mo* A“ A^) ■«" *• m“ * * /\ 0-) m 

where 7((m)-*- 1 if m-s*- 00 . When and ^ are fixed the results de- 
pends only on the product Nq • A*A eff .The eq. / 7 / describes the 
asymptotic behavior of expression / 5 / shown as solid curves 
on fig. 1 .By definition the effective atomic number for a mixed 
composition is: 

<=■ y/s . £- 


A, 


which we can evaluate comparing / 7 / with experimental data. 
For _ tf = 1.7, 5 =0.73,Ng = 0.il3±0.12 J A eff = 0.15i0.02 (of them onl 
N n is model dependent ('5.6>lthe new value nf ns = 


If 


i *- s 

XT 


A: , 

lA; .. 


78/ 


0 is model dependent (5,6))the new value of A eff is A e ff= 
3. 9*1. 5 (using experimental value of I(m) at m=10). 


only 


Table 1 


5* (Gev/n) 

A 

1 

4 

9 

14 

28 

56 

9^ 

.939 

.055 

.0009 

.0035 

.0011 

.0003 

CC1 


1.7 

1.7 

1.7 

1.7 

1.7 

1.7 

EC 2 


1.7 

1.7 

1.7 

1.7 

1.6 

1.6 


Table 1 shows two suggested examples of chemical composi- 
tion, 9‘t -fraction of nuclei Ai with a given energy per nucle- 
on (at 1 Gev/n) .Number spectra for these two cases calculated 
using eqs. 71-5/ shown in fig.l by dashed lines. A good fit to 
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experimental data is provided by CC1 which corresponds to the 
case of a constant chemical composition (a standard one) in 
all energy range from 1 Gev/n.The A e ff for CC1 according to 
/ 8 / is A ef £=4 . 8 , which is about 20% bigger than obtained using 
■'No-approximation" at m=10;this is because f(Eg) did not 
reach asymptotic behavior at m=10 for heavy primaries .Diffe- 
rent multiplicities m represent different primary energies 
but also different selection of A. Table 2 shows , assuming CC1 
and eqs . / 1 — 5 / , relative contributions of different A, also 
their mean energy,to different m,q is relative contribution. 


Table 2 



m = 

1 

2 

3 

5 

10 

30 

A = 1 

E = 

6.1 

25 

86 

240 

520 

2500 


q = 

.83 

.64 

.47 

.32 

.22 

. 16 

A = 4 

E = 

20 

44 

94 

170 

340 

1600 E in 


q = 

.14 

.26 

.33 

.33 

.26 

. 2 Tev/nucleus 

A = l4 

E = 

53 

78 

106 

160 

250 

1050 


q = 

.02 

.07 

.13 

.2 

• 23 

.2 

A = 56 

E = 

170 

200 

210 

220 

245 

730 


q = 

.005 

.03 

.07 

.15 

.3 

.44 


4 . Conclusions . The Baksan muon number spectrum at E jJ >E t j ;i = 
=220 Gev is well explained suggesting a constant chemical 
composition (the same as at lGev/n) till'~20 Tev/n. 

The integral exponent of the heavy primaries power law 
spectrum should be ^ = 1.7 1 0.1 till 10l5 ev/nucleus. 
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